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An exact solution to the Friedmann equations with a string inspired phantom scalar 
matter field is constructed and the absence of the "Big Rip" singularity is shown 
explicitly. The notable features of the concerned model are a ghost sign of the kinetic 
in ■ term and a special polynomial form of the effective tachyon potential. The constructed 

solution is stable with respect to small fluctuations of the initial conditions and special 



deviations of the form of the potential. 



(N 

^ : 1 INTRODUCTION 
o 

^ ■ The combined analysis of the type la supernovae, galaxy clusters measurements and WMAP 
^! (Wilkinson Microwave Anisotropy Probe) data brings out clearly an evidence of the accel- 
O I erated expansion of the Universe The cosmological acceleration strongly indicates 

I that the present day Universe is dominated by a smoothly distributed slowly varying cosmic 
■ fluid with a negative pressure, the so-called dark energy |S]-|H] (alternative approaches are 
presented, for example, in 0). 
^ , To specify different types of cosmic fluids one usually uses a phenomenological relation 

H I between the pressure density p and the energy density g, corresponding to each component 
of fluid 

p = wg, 

where w is the equation-of-state parameter or, for short, the state parameter. A component 
with negative w corresponds to the dark energy. Contemporary experiments, including 
WMAP, give strong support that at present time the dark energy state parameter is close 
to —1 jl], [ini-[lll- In particular, it follows from the current observational bounds j2] that 
w belongs to the interval: 

-1.46 <w < -0.78. 

at the 95% confidence level. 
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From the theoretical point of view the above-mentioned domain of w covers three essen- 
tially different cases: w > w = —1 and w < ~1. 

• The first case, w > —1, is achieved in quintessence models PS]-|221, which are cosmo- 
logical models with a scalar field. Such types of models are quite acceptable, but there 
is a question of an origin of the scalar field. To comply with astronomical experimental 
data this scalar field should be extra light and hence it does not belong to the Standard 
Model set of fields |2S]- 

• The second case, w = —1, is realized by means of the cosmological constant. This 
scenario is admissible from a general point of view except for a problem of an order 
of the magnitude of the cosmological constant. It should be 10^^*^ times less than the 
natural theoretical prediction jH]. 

• The third case, w < —1, is called a "phantom" one and can be realized due to a scalar 
field with a ghost (phantom) kinetic term. In this case all natural energy conditions are 
violated and there are problems of instability at classical and quantum levels [^-|44j. 
Since experimental data do not contradict with a possibility of w < — 1 and moreover 
a direct search strategy to test the inequality w < —1 has been proposed it is 
interesting to find a noncontradictory model with the condition w < —1 satisfied^. 

Note the K-essence models ^H] can have both w < —1 and w ^ — 1. At the same time a 
dynamical transition from the region w ^ — 1 to the region w < —1 or vice versa is forbidden 
under general assumptions [Hj and is possible only under special conditions |48j. 

Let us recall that in models with a constant state parameter w, less than —1, and the 
spatially flat Friedmann metric the scale factor tends to infinity and therefore the Universe 
breaks down at a finite moment of time. This problem is known as the "Big Rip" [2^1; see 
also [HZj. The simplest way to overcome this difficulty in models with w < — 1 is to consider 
a scalar field with a negative time component in the kinetic term f2ll2H|- The appearance 
of a such summand is possible due to quantum effects [IHIEII- However, all proposed models 
insist of the instability problem. 

A possible way to evade the instability problem for models with < — 1 is to yield 
a phantom model as an effective one, arising from a more fundamental theory without a 
negative kinetic term. In particular, if we consider a model with higher derivatives such as 
0e~'~'0, then in the simplest approximation: 0e~^0 ~ 0^ — 000, such a model gives a kinetic 
term with a wrong (ghost) sign. It turns out, that such a possibility does appear in the string 
field theory framework jH] . Since the concerned model is a string theory approximation, all 
stability problems related to a model with higher derivatives are discarded. 

Our goal in the present paper is to construct an analytic solution in a polynomial model, 
which is close in some approximation to a model arising in the string theory, namely in 
the theory of a fermionic NSR string with GSO— sector. The scalar field is an open 
string theory tachyon, which, according to the Sen's conjecture [21], describes brane decay, 
at which a slow transition to the stable vacuum, correlating with states of the closed string, 

^Let's note the model in which w < — 1 is a result of axion consideration. 
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takes place. Other dark energy models, which use the brane-world scenarios, are presented 
in ISO]. 

The notable feature of the model is a ghost sign of the kinetic term and a special form 
of the potential: 

vi<f^) = l{i-<p'Y + j^f{^-<p'Y. (1) 

In the first approximation only the first summand in the right hand side of (0) can 
be obtained from the open string field theory in a flat metric (which corresponds to 
nip = oo). The second summand allows us to save the form of the interpolating solution 
obtained in the fiat metric for the case of the Friedmann Universe, i.e. at an arbitrary 
rrip. Note, that the appearance of higher powers of is possible if we take into account 
the heavier perturbations of the string (52j, but we can not expect to obtain the explicit 
values of the coefficients of the potential (QJ from the string theory. The performed in this 
paper stability analysis of solutions points out that the constructed solution is an attractive 
one and there are no essential changes in the behavior of the solutions for some variations 
of the potential parameters. This gives reasons to assume, that the qualitative behavior of 
the physical parameters H,w,q in the model with the potential (Q) under certain conditions 
becomes close to the behavior of the corresponding ones in the model with a string theory 
potential. 

The constructed solution has no singular points and becomes the de-Sitter solution at 
large time. It does not contradict to the general theorems of the relativity theory because 
the energy conditions are violated in the considered model (this question is discussed in ^53^). 

The paper is organized as follows. In Section 2 we describe the cosmological model, 
which is an approximation for the string inspired model, and write down the solutions to 
the Friedmann equations with the special polynomial potential. In Section 3 we consider 
the phantom dynamics in detail and find out that our solution describes an accelerating 
Universe. In Section 4 we study the stability of the solution with respect to variations of 
initial data and peculiar variations of a form of the potential. In Section 5 we do sum up 
the obtained results and point out the ways for generalization of the considered model. 



2 MODEL 

2.1 Equations of motion 

Let us consider the gravitational model with a phantom scalar field and the spatially flat 
Friedmann metric. Since the origin of the phantom field, considering in this paper, is con- 
nected with the string field theory we include the typical string mass Mg and a dimensionless 
open string coupling constant Qo in the action. The action is: 

S = I d'x^g (^^i? + 1 [+\g^^d,<pdA - V{<P)^ ) , (2) 

where Mp is the Planck mass, g is a. spatially fiat Friedmann metric: 

ds'^ = — dt^ + a^(t) {dx\ + dx\ + dx^ , 
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where a{t) is the scale factor. The coordinates {t, Xi) and the field (p are dimensionless. 
Hereafter we use the dimensionless parameter rrip = g'^M^ / . 
Independent equations of motion are as follows 

3H^ = ^g, + 2H = -\p. (3) 

Here H is the Hubble parameter: H = a{t)/a{t) and the dot denotes the time derivative. If 
the scalar field depends only on time, i.e. = then 

Q= -\^' + V{<P), P=-^0'-m- (4) 
We recast system (jH)) to the following form 

1 :o If 1 : 



Besides of equations (jSJ one can obtain from the action ^ an equation of motion for the 
field 

+ 3i/0-l^; = O, (6) 

where = This equation is in fact a consequence of system Q. Expressing H through 
(j) and 0, we obtain 



+ sign{H)^U-h^ + V{4>) - = 0. (7) 

Here sign(i7) factor selects the proper branch of the square root. The expanding Universe 
corresponds to H > 



2.2 Connection with the String Theory 

At present time one of the possible scenarios of the Universe evolution considers the Universe 
to be a D3-brane (3 spatial and one time variable) embedded in higher-dimensional space- 
time. D-branes do naturally emerge in the open string theory (in our case we consider 
fermionic open strings containing both GSO+ and GSO— sectors). The D-brane in question 
is unstable and does evolve to the stable state. This process is described by the open string 
dynamics, which ends are attached to the brane (see reviews [^2] and references therein). 
If only the lowest excitation — the tachyon — is taken into account then the D-brane 
dynamics are described by an action of the open string tachyon. There are two common 
ways to describe the tachyon behavior: Dirac-Born-Infeld (DBI) approximation and the 
level truncation scheme in the covariant string field theory [SHISSlEni- The level truncation 
method advocated itself checking the Sen's conjecture and, therefore, it is natural to use it to 
analyze the dynamics of the D-brane, see |F71 158| 1^. The action for the tachyon resulting 
from the fermionic string field theory in the approximation of the slow-varying auxiliary 
field has the form 
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where Lp is the dimensionless tachyon field and the coordinates and the constant go are also 

dimensionless. The tilde means an action of the non-local operator exp ^— log ^^^j 

on the field y?. The subscript "fiat" hereafter indicates the fiat background metric is used. 
Note that in the bosonic string case the interaction is cubic. For the space homogenous 
configurations only the dependence on time r is left and the action takes the form 





2 \dT ) 256 

Introducing the notation = j:^exp ^log ^y^^j^ ^ and t = r/2, we rewrite the equa- 
tions of motion obtained from the last action in the following form 



a3 
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where d = "i^oS^^^^ > 0- f^^t this equation is integral one since the exponent 

containing an infinite number of derivatives can be replaced with an integral operator jHZj 
in some class of fields as follows 

Note, that equation (jSI) without the first summand is an equation for the p-adic string (in our 
case p = 3) jHO]- In the latter case the existence of a solution interpolating between vacua 
has been proved jHJ. Previously this solution has been numerically constructed in [U^ . 

Equation (jH)) was investigated in [SHI EHl and an existence of a rolling solution, i.e. a 
solution interpolating between the two vacua of the potential 

VU<P) = \{l~<P')' (9) 

was shown numerically. It is remarkable that the late time behavior of the rolling solution 
can be effectively described by a lagrangian with a ghost sign of the kinetic term and the 
same potential. The corresponding action has the form 



^Sflat= J dti^-^<p^--{i-<pr)- (10) 

As known, the equation of motion derived from the action ()10|): 

+ 20(l-02)=o (11) 
has an interpolating solution — a kink 

(f){t) = tanh(t). (12) 
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2.3 Construction of solution 



System of equations © with some polynomial potential V{(j)) is not integrable. The nu- 
meric analysis of the Friedmann equations with the potential has been performed in the 
paper [Hj. It is known ( see, for example [33 IHSl) that it is always possible to find the 
potential V{(f)), if either or H(t) are given explicitly. 

It is worth to note that our model has common features with models considered in brane- 
world scenarios (see Appendix). Following [63^, we assume, that H(t) is a function of (pit), 
called superpotential, that is H(t) = W{(f)(t)). Using equality H = W^cj), we obtain from 
system dHJ: 

^ = 2mjiy;, (13) 
V = 2mjl^;' + 3mlW^ (14) 

Let us find a potential, which corresponds to the field (p of type (|T^ . This function 
satisfies the following equation 

= 1-0', (15) 

hence, from ()13|) it follows 

where C is an integration constant. Once W{(j)) is known we obtain the potential V{(f)) from 
()14j) . Different values of C correspond to different V{(f)). The requirement that V{(j)) is an 
even function is equivalent to the condition C = 0. In this case 



n0)=^(i-0')VY^0'(3-0'r- (16) 



p 



It can be straightforwardly verified that if (j){t) is a solution to (jl3p . (jl4j) with some 
potential V{(j)), then 4>{—t) is a solution as well. For example, 02(t) = tanh(— t) = — tanh(t) 
also is a solution to (fT^ . (fT^ with the potential (fT^ . Indeed these two solutions correspond 
to the one and the same function H(t), but to the two different functions W{(j)). 

We have constructed the potential V, using an explicit form of a solution. Now we 
consider more general problem and look for solutions, starting with the the following two 
requirements: (j){t) is a smooth function with nonzero asymptotics 0(±oo) = ±B and the 
superpotential W{(j)) is a polynomial in (p not more than the third degree. In the first 
condition we can put B = 1 without loss of generality. The second condition guarantees 
that the potential V{(p) will be a polynomial in (p and its degree is not more than six. The 
above conditions allow us to rewrite equation ()13|) in the following form: 

= a + /?0 + 702. (17) 

From the asymptotic conditions it follows 

7 = -a, p = 0. (18) 



6 



Taking into account (fTHj) we solve equation (fT7j) and obtain 

=tanh(a(t-to))- (19) 

We do not consider another solution: (pit) = cotli(a(t — to)), because it is not a smooth 
function. Thus, up to a scaling and a time shift we obtain a solution, which coincides with 
the solution in the flat metric. It is evident that the potential obtained from (fT^ coincides 
with ()16p up to an overall factor. It is interesting that under above-mentioned conditions 
on a solution and a superpotential, the sixth degree potential is a minimal possibility for 
existence of an interpolating solution. One can also see, that any solution to the equations of 
motion in flat space with the fourth degree potential is a solution to the Friedmann equations 
either with the sixth degree polynomial potential or with a nonpolynomial potential 64 . 

It is interesting to note that we have built an exact solution without any approximation. 
One of the standard approaches to the analysis of cosmological solutions in the presence of 
a scalar field is the so called slow-roll approximation which technically means a neglecting 
of the terms in equation (jHI) and 0^ in the expression for the energy (@]). However, from 
(fT3j) it follows that for our solution the following relations take place: 

= -20(l-0'), 0^= (1-0^)'. 

It is evident that the slow-roll approximation in this case is valid only when is close 
to 1. The above-mentioned approximation describes our solution at very large times and 
can not be used for its description in the beginning (close to zero) moments of time. 

The Hubble parameter for the constructed solution is given by 

This parameter goes asymptotically to l/(3mp) when t — * oo, or equivalently 0^1. Once 
H{t) is known one readily obtains the scale factor 

e^'/(^2m2) , , /cosh(t)2-l\ ^ ^ 

a = an , = ao(cosh(t))'^'> exp , , ,^ , (20) 

"(1 _02)i/(6m2) w; ^ Vl2m2cosh(t)2y ^ ^ 

where ao is an arbitrary constant. 



3 COSMOLOGICAL CONSEQUENCES 

3.1 Acceleration 

The function a{t) has the following asymptotic behavior at large t: 

lim ait) ~ e'^'^v . (21) 

t— »oo 

It follows from formulae (|2Up and (j2H) that the Universe expands with an acceleration. The 
deceleration parameter is negative and equal to 

aa 2 5 + cosh(4t) 



<l{t) = - 



d2 sinh(t)2 2(2 + cosh(2t))2' 
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The corresponding plots are drawn in Fig. ^ (Hereafter we assume rrip = 1 for all plots if 
something else is not specified explicitly). The "jerk" parameter j = 'do? /a^ and the "kerk" 




Figure 2: "jerk" parameter (left), "kerk" parameter (center) and its fine structure (right). 

Substituting the obtained solution (fT^ and the potential (fT^ into expressions for the 
pressure and energy densities, we obtain 



= - (1 - 0^)' - (3 -<P')\ m = (3 - 



12mp ' 12mH 

Therefore, the parameter of state w is given by 



p 9 (1 - (p^ f 9 cosh(t)^ 

w = - = -l- 12m! ^ = -1 - 12m? ^ ^ 



Q ' " ^02(3_02)2 ^sinh(t)2(l + 2cosh(t)2)2' 

Note that the function w{(f)) is less than —1 for |0| 7^ 1 and is equal to —1 for |0| = 1. 
Point 0=1 corresponds to infinite future. The plots for the pressure density p{t), the energy 
density g(t) and the parameter of state w are presented in Fig. |21 

Thus, in the considered model of the accelerated expanding Universe the phantom field 
describes the dark energy with the state parameter w < —1. 
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Figure 3: p(t) (left), g(t) (center) and w(t) (right) 



3.2 Evolution of the exact solution to equation ([70 

Solution (fT^ to equation ((Tj) is a function describing an interpolation of the field from 
point —1 to point 1 (and not viceversa) with zero initial and final velocities in the potential 



(fT(3|) with a friction proportional to sign(if) y — + V{(f)), which depends on coordinate 
and velocity. Factor sign(if) makes the friction negative for negative 0, i.e. the particle is ac- 



/{2V3: 



rrir. 



celerated by it. Indeed, the expression for friction coefficient is equal to (3 — 
on the solution for both positive and negative 0. 

Let us discuss on the evolution of the exactly constructed solution. A phantom field 
evolution is equivalent to an evolution of the ordinary field in the fiipped potential. The shape 
of potential V essentially depends on value of nip, see Fig. HI where the flipped potential —V 
is plotted for different values of nip. 

To analyze the extrema, let us note the following expressions 



V{0) 



1 

2' 
3 



V{±1) 



3m2' 



(1-0^) ((3-4mJ)-0^) 



Arrip 



2 {2ml - 1) 



(22) 



It is evident from 



that extrema 



±1 and 



are independent of nip, whereas 



other extrema are functions of nip. Moreover, a particular value of nip determines whether 
a specific extremum is a maximum or a minimum. Possible values of nip can be separated 
in the following domains which determine a structure of the fiipped potential —V extrema. 



• For rup = 00 V{(j)) is a well known double well potential. 

• For I ^ < 00 the fiipped potential —V, being a 6th degree potential is similar 
graphically to a fiipped double well potential: it has a local minimum in the point 
= and maxima in the points = ±1. 

• For I < nip < I point = in the fiipped potential becomes a local maximum two 
more minima appear in the interval (—1,1), developing two wells and one hill on the 
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Figure 4: Potential —V for different values of nip. The following nip correspond to the 
curves (top-down direction): = oo, = 1, = 1/2 and = 1/4 in the left plot and 
nip = 17/24 and = 7/12 in the right plot. Two vertical lines = ±1 are used as a grid. 



way of our field during an interpolation from point —1 to point 1, see the right plot 
in Fig. ^ In more details: 

— For I ^ nip < I points = ±1 are not below the point = in the flipped 
potential. 

— For i < nip < I points = ±1 are below the point = in the flipped 
potential. This means that our field (p starts at the point (p = —1 with zero initial 
velocity, rolls down to the well and then climbs up to the hill which is above its 
initial location. This is counterintuitive while the factor sign (if) in the friction 
coefficient is not taken into account. The latter means the friction is negative 
for negative 0. It seems one can disprove the last statement since it is natural 
to inverse the time and consider the backward motion. It is easy to resolve this 
contradiction because in the case of an inverse time function H would have an 
opposite sign. 

• For nip = i the flipped potential has the only extremum which is maximum in the 
origin and points = ±1 are inflection points. 

• For TUp < I two minima, coming from the last multiplier of the expression V fl22|) . 
go beyond the interval — 1 < < 1 and become maxima. Points = ±1 become 
minima of the flipped potential. Point = remains the maximum. In this case field 
starting at point = — 1 climbs up to the hill from the very beginning and then rolls 
down and stops at the point = 1. Such a behavior seems incredible and becomes 
real due to a negative friction for negative 0. 
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Investigating the cosmological evolution of the phantom field we use only the a part of 
solution (fT^ starting at the time t = tj„ > 0. At this time the field is at the point (pin = 
^{tin) = tanh(tj„) and one has to supply an initial velocity 0(ti„) = 1 — 0f„ = 1/ cosh^(tj„). 
For large (weak coupling to the gravity) the phantom field uses an initial kinetic energy 
to climb up to the hill at = 1. Increasing a coupling to gravity, i.e. decreasing we lower 
the height of the hill at the point = 1 of the fiipped potential. For rrip = 2/3 heights of 
the hills at points = and 0=1 become equal and all initial kinetic energy is spent for 
the work against the friction. For < 2/3 the friction becomes stronger and the particle 
has the kinetic energy which is exactly enough to reach the point = 1 which is now lower 
than the hill at point = 0. 



4 STABILITY OF SOLUTIONS 

4.1 Variation of initial data 

Let us consider the behavior of the solution to system in the neighborhood of our exact 
solution 

0o(t) = tanhm, Ho(t) = ^tanh(t) ( 1 - -tanh^^ 

2m^ \ 3 

with the aim to analyze its stability. Substituting 0(t) = 0o(t) +e0i(t) and H{t) = Hoit) + 
eHi{t) in system Q, we obtain in the first order of e the following equations: 

Hi = ^fl -tanh(t)2)0i, 

(18 - 24ml + - 1) tanh(t)2 + 6 tanh(t)^) tanh(t) 



^ ~ 12m2 (1 - tanh(t)2) 

(3 -tanh(t)2)tanh(t)^ 
1 - tanh(t)2 ^' 

System (PHj) has the following solution: 

01 (t) =2mJCi (l-tanh(t)2) + 



(23) 



, 2-- 



2„ 2J{t) + (cosh(2t) - l)(cosh(t)) ™^ eV2'"^(-=h(2t)+i)^ 
+ 2^p^2 ^^.U9,^M^ ' (24) 



Hi{t) = Ci (1 -tanh(t)2)^ - 



cosh(2t) + 1 
AmlC2J{t) 
;cosh(2t) + l)2^ 



where 



ft 2 1 

J(t) = / sinh (r) (cosh(r))^~^/"'' (2 (2mJ - l) cosh(r)2 - l) e^'-^-'^f^'" rfr. 

If C2 = 0, then Hiif) and 0i(t) are bounded for all t G (—00,00) and for all values of 
nip. In this case the solution can be presented in the following form: 

0(t) = 0o((l + eCi)t) + 0(£2), Hit) = Hoiil + eCi)t) + 0{e'). 
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Let us consider the case C2 7^ 0. It is easy to see that if > 1/2 then 0i(t) goes to 



infinity as t 
(01) that 



00 and, therefore, our solution is not stable. At ml = 1/2 we obtain from 



Hi{t) = (tanh(t)2 - 1)\C, - C2J2 

Mt) = 



(tanh(t)2 - 1) (Ci - C2J2) - -Cae-^'^^'^W 



where J2 = /oe"*^°''W'/2^anh(r)rfr. Thus, and (pi are bounded functions. For < 1/2, 
01 (t) and -ffi(t) are bounded functions as well. Thus we obtain that our solutions are stable, 
the first corrections are bounded, if ^1/2. We remind, that for < 1/2 points = ±1 
are minima, and for = 1/2 these points are infiection ones. 

To study a behavior of solutions with initial conditions much distinct from the initial 
conditions of the exact solutions and to answer the question whether the obtained trajectory 
is an attractive or repulsive one, we construct phase portraits for various values of m^. 

Note, that H{t) is a real function. It gives a restriction on the maximal value of 4>{ty, 
namely, for all the inequality that 0^ ^ 2V^(0) is hold. For example, let 0(0) = then only 
solutions with 0(0) ^ 1 have the physical sense. In other words our exact solution starts 
from this point with the maximal possible speed. 

We can not find exact solutions to system © with initial conditions 0(0) = and 
10(0)1 < 1, and thus we present the results of numeric calculations. 

In Fig. El phase portraits for rrip > 1/2 are presented, i.e. in the domain of instability 
of our solution. One can see on the graphs that evolution of solutions comes to the end at 
points of a minimum of fiip potential —V, namely at zero for = 100 and = 1 and at 
the point = ^/2/3 ^ 0.816 for = 7/12 (see. Fig. ij). 
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Figure 5: Phase portraits for different values of > 1/2. Thick lines strokes depict solutions 
with the initial condition 0(0) = 1/2. 

Let's consider phase portraits at ^1/2, presented in Fig. IHl In the case = 1/2 
when = ±1 are inflection points, a solution with the initial speed 0(0) = 1/2, as well as 
the exact solution tends to the point 0=1, however coincides with it only in this point. For 
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m-p = 1/4 and = 1/100 we see, that the numeric solution is close to the exact solution 
for (f) < 1, and the less the less is the value of 0, at which solutions with initial speeds 



0(0) = 1/2 and 0(0) = 1 get an identical speed. 



m 
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Figure 6: The phase portraits at ^ 1/2. Thick lines strokes depict solutions with the 
initial conditions 0(0) = 1/2 and 0(0) = 1. 

Thus both behavior of the first corrections and phase portraits show that the obtained 
exact solution is stable for ^1/2. 

4.2 Variation of the form of the potential 

In this subsection we investigate an influence of variation of the coefficients in the second 
summand of the potential ()16|) on the behavior of solutions to system (0). Consider a 
dynamics of system (jSJ with the potential 

Parameter y is introduced so that the points = ±1 remain the extrema of the potential. 
Note, that up to a constant shift and an overall factor multiplication this is the only possible 
modification of an even potential (|T6|) . which leaves the points = ±1 to be extrema. 
Equation (0) differs from the kink equation of motion (fTT|l by the friction like term and an 
extra proportional to l/rUp term in the potential. Parameter y enables us to separate the 
effects of the friction and an extra term in the potential and to answer the question which 
one provides the stability of solution. 

To analyze stability it is convenient to plot phase portraits for different values of pa- 
rameters rup and y and to render them in a table as it is done in Figs. U\ and |H1 Guided 
by the plotted phase portraits and an analysis performed in Sections 3 and 4, we conclude 
the stabilization of solutions is due to an extra term in the potential. Namely, the obtained 
solution is attractive one for m^/y ^1/2 irrespectively the numeric factor in the friction 
term. 
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Figure 7: Phase portraits for rrip ^1/2 and various y. 
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5 Conclusion 



In the present paper we have performed an analysis of an exactly solvable model of acceler- 
ating expanded Universe dominated by the dark energy. The state parameter of the model 
in the consideration is always less than —1 and tends to —1 when t —>■ oo. Such a behavior 
of the state parameter leads to an absence of the Big Rip singularity in our model. This kind 
of singularity does exist in the model with a constant w < —1. Note, that these properties of 
the state parameter are kept after including in our model an interaction with the cold dark 
matter [65J. It is shown also in this paper that our solution is stable with respect to small 
fluctuations of initial data when ^1/2 and small fluctuations of the form of the poten- 
tial. The simplest way to yield an exactly solvable model with w crossing the cosmological 
constant barrier w = —1 consists in an inclusion of an additional scalar fleld [66]. 

Authors are grateful to A.E. Pukhov and I.V. Volovich for useful discussions. The work 
is supported in part by RFBR grant 05-01-00758, I.Ya. and A.K. are supported in part by 
INTAS grant 03-51-6346 and Russian President's grant NSh-2052.2003.1, the work of S.V. 
is flnancially supported in part by Russian President's grant NSh-1685.2003.2 and grant of 
the scientiflc program "Universities of Russia" 02.02.503. 



Appendix. Cosmo-phantom versus 5-dimensional grav- 
ity models 

Our model as whole and system in particular are in close analogy with models consid- 
ered in brane-world scenarios. Let us compare system (0) with the system of the Einstein 
equations in the model of gravity interacting with a single scalar fleld in flve-dimensional 
space-time [^3 ISZj • The action is 

S = J d'xdr ^\det - ^r''d,4>du4> - ^(0) j , 

where ^^jy is deflned as follows 

ds^ = e^^^''\-dxl + dx\ + dxl + dxl) + dr"^. 
Assuming the scalar fleld depends only on r, the independent equations of motion are 

H' = (25) 

3^^ = -y (26) 

where H = 4A'/3 = |^ and 0' = These equations correspond to (0) with = —9/16. 
Surely, such a choice of the coupling is not physical, but we see that the method described 
in jnSl to transform (0) into the flrst order differential equation works in our case as well, i.e. 
mathematical method, which simplifles (|25|) and (|26p. does not depend on a value or sign 
of rni. 
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